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The Unimodular Intersection Problem 
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Abstract 

We show that finding minimally intersecting n paths from s to f in a 
directed graph or n perfect matchings in a bipartite graph can be done in 
polynomial time. This holds more generally for unimodular set systems. 
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1 Introduction 

Let G be a directed graph on d edges with two distinct vertices s and t and let n 
be a positive integer. Given n paths from s to t, an edge is critical if it is used by 
all paths. We wish to find n paths from s to t with minimum number of critical 
edges. One motivation for this problem, adapted from mi, is as follows. We 
need to make a sensitive shipment from s to t. In the planning stage, n paths are 
chosen and prepared. Then, just prior to the actual shipment, one of these paths is 
randomly chosen and used. An adversary, trying to harm the shipment and aware of 
the prepared paths but not of the actual path that is finally chosen, will try to harm 
a critical edge, used by all paths. Therefore, we protect each critical edge with high 
cost, and so our goal is to choose paths with minimum number of critical edges. 

Here we study and solve a more difficult problem, described as follows. Following 
PQ, for 1 < fc < n, call an edge k-vulnerable if it used by at least k of the n paths. In 
[3] it was shown that finding n paths with minimum number of 2-vulnerable edges 
is NP-hard, and in [I] it was shown that for fixed n and k, finding n paths with 
minimum number of ^-vulnerable edges can be done in polynomial time. 

Here, we assume that n is variable, and want to find n paths with lexicographically 
minimal intersection, that is, which first of all minimize the number of n -vulnerable 
edges, then of (n — Invulnerable edges, and so on. More precisely, given n paths, 
define their vulnerability vector to be / = (fi, ■ ■ ■, f n ) with f k the number of k- 
vulnerable edges. Then n paths with vector / are better then n paths with vector g 
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if the last nonzero entry of g — f is positive, denoted / -< g. In particular, n paths 
with lexicographically minimal vector have a minimum number of critical edges. 

Let A be the vertex-edge incidence matrix of G and let b be the vector in the 
vertex space with b s = —1, b t = 1, and b v = 0 for any other vertex v. For x G M. d let 
M := Yli=i \ x i\- For • • •, x n G R d let x 1 A • • • A x n G R d and x l V • • • V x n G be 
the entry-wise minimum and maximum of x 1 ,..., x n , and dehne /(x 1 ,..., x n ) by 

/^(x 1 ,..., x n ) := | V {x n A • • • A x lk : 1 < i\ < ■ ■ • < < n}\ , k — 1,..., n . 

So for x 1 ,..., x n G {0, l} d we have that fkix 1 ,..., x n ) is the number of k-vulnerable 
components 1 < i < d, in which at least k of the vectors x 1 ,... ,x n have a one-entry. 
Then our problem is the following nonlinear integer program defined by A, 6, n: 

lcxmin{/(x 1 ,..., x n ) : x k G {0, l} d , Ax k — b , k — 1,..., n} . (1) 

In this article we solve this problem, even in the more general context where A is 
an arbitrary totally unimodular matrix. We establish the following theorem. 

Theorem 1.1 The lexicographic problem ijj\) over any totally unimodular matrix A, 
any integer right-hand side b, and any positive integer n, is polynomial time solvable. 

Another consequence of this theorem is the following. Let G be a bipartite graph 
on d edges and let n be a positive integer. Given n perfect matchings in G, an edge 
is vulnerable if it is used by at least k matchings. The problem of finding n perfect 
matchings with lexicographically minimal intersection is again a special case of our 
theorem with A the vertex-edge incidence matrix of G and b the vector in the vertex 
space with b v — 1 for every vertex v. So it can also be solved in polynomial time. 

In the next section we prove our theorem, proving on the way a core lemma 
which is of interest in its right. In the last section we discuss possibilities of some 
further generalizations, solve the simpler problem of minimizing f n (x 1 ,...,x n ) = 
lx 1 A • • • A x n \ only more efficiently, and describe some NP-complete variants. 

2 Proof 

We begin with an algorithmic version of a decomposition result of [0]. 

Lemma 2.1 There is a polynomial time algorithm that given any positive integer n, 
any totally unimodular matrix A, any integer vector b, and any vector 
x G {0,1,..., 7i} d satisfying Ax = nb, finds vectors x 1 ,... ,x n G {0, l} d satisfying 

n 

x k = x , Ax k = b , k = 1 ,... ,7i . 

k =1 
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Proof. We use induction on n. For n — 1 take x 1 = x. Next consider n > 1 and 
the system 0 < < Uj < \~ x j\ < 1, j = 1,..., d, Ay = b. Since -x is a real 

solution and A is totally unimodular, there is also an integer solution x n G {0, l} d , 
which can be found in polynomial time by linear programming. Set z := x — x n . 
Then Az = (n — 1 )b, and for all j, if Xj = 0 then xf = 0 and if Xj = n then 
xf — 1, so z G {0,1,..., n — l} rf . Therefore, by induction, we can find in polyno¬ 
mial time vectors x 1 ,... ,x n ~ 1 G {0, l} d satisfying Ylk = = z and Ax fc = b for 
all k. Then x k = z+x n = x and the induction follows, completing the proof. □ 


Next, we make the following definition, which is central to what follows. The 
compression of x — (x 1 ,..., x n ) G {0, l} dn is x = (x 1 ,..., x n ) G {0, l} dn dehned by 

x k = V{x n A • • • A x lh : 1 < i\ < ■ ■ ■ < ik < n} , k — l,...,n. 

Thus, the vulnerability vector of x — (x 1 ,..., x n ) G {0, l} dn is f(x) = (lx 1 ),..., |x n |). 

Call c = (c 1 ,..., c n ) G M. dn nondecreasing if c 1 < ■ ■ ■ < c n . For c, x G M dri put 
cx := J^! =1 c k x k . We have the following lemma which is of interest in its own right. 

Lemma 2.2 There is a polynomial time algorithm that given any positive integer 
n, any totally unimodular matrix A, any integer vector b, and any nondecreasing 
c = (c 1 ,..., c n ) G 7L dn , solves the following nonlinear integer programming problem, 

minjcx : Ax k = b , k — 1,... ,n, x = (x 1 ,..., x n ) G {0, l} dn } . (2) 


Note that the objective cx is nonlinear since c acts on the compression of x. 
Proof. First, we need a couple of properties of the compression. Note that it satisfies 

1 if k < £"= i x \ 


x k = 


0 otherwise 
Therefore, for all x,y,z G {0, l } dn , we have 

n n 

£** = £**, 


, k — 1,..., n , i — 1,..., d . 


k =1 


k =1 


and , if , then y = z . 

k=l k =1 

Now, consider the linear integer program 

n 

min{cx : A ^ x k = nb , x = (x 1 ,..., x n ) G {0, l} dn } . 


( 3 ) 

( 4 ) 


( 5 ) 


k =1 


As [A,..., A] is totally unimodular since A is, and b is integer, we can solve program 
(l5l) in polynomial time using linear programming. If it is infeasible then so is program 
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© and we are done. So assume we obtain an optimal solution y G {0, l} dn to ([5]). 
By Lemma [2.11 applied to Y^k=\ V k we can find in polynomial time z G {0, l} dn with 

n n 

Az k = b, k = l,...,n. 

k =1 k =1 

Clearly z is feasible in 02]) • We claim it is also optimal in 02]) • To prove this, we 
consider any x which is feasible in 02]), and prove that the following inequality holds, 

cz = cy < cy < cx . 

Indeed, the first equality follows since z k = J2k=i V k an d therefore z = y 

by 03) • The middle inequality follows since c is nondecreasing. The last inequality 
follows since Ylk=i xk = J2k=i xk by © and therefore x is feasible in ©. □ 

We now proceed with the proof of Theorem II.It Define a nondecreasing c G Z dn , 
c k := (d+ l) fc_1 , k — l,...,n, i — 

Consider any x, y G {0, l} dn , with vulnerability vectors / = (lx 1 ),..., |x n |) and 
g = fly 1 !,..., \y n \). Suppose / is lexicographically smaller than g, that is, / -< g. 
Let r be the largest index such that \x r \ ^ \y r \. Then \y r \ > |x r | + 1. We then have 

n 

cy-cx = ^(d + I) 1 * -1 (\y k \ - \x k \) 

k=1 

> ^(d+if-qiyi-^D + ^+ir 1 

k<r 

> (d + iy- 1 -J2 d ( d + i) fc_1 > o. 

k<r 

Thus, an optimal solution x 1 ,... ,x n for © is also optimal for our original prob¬ 
lem ©, and can be found by Lemma 1X21 in polynomial time . Theorem II.ll follows. □ 


3 Remarks 

Gijswijt [2j established a decomposition result like the one we recalled in Lemma [XT] 
for the more general case of A being only nearly totally unimodular , i.e., A = A'+ca T 
with the matrix arising from adding row a T to A' being totally unimodular and the 
vector c being integer. Using this result, one can, for instance, see that by the 
algorithm described in Section [2] the lexicographically minimal intersection problem 
cannot only be solved in polynomial time for structures like s-t paths or matchings 
in bipartite graphs, but also, for example, for stable sets in circular arc graphs. 
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Next we discuss the simpler problem of minimizing ■ ■ ■, x n ) = (a; 1 A- • • Ax n \ 

only, which in the paths problem is the number of critical edges, that is, the problem 

min {Jo; 1 A ■ • • A x n \ : x k G {0, l} rf , Ax k = b , k = 1,..., n} . (6) 

It can be solved in polynomial time as a corollary of Theorem 11.11 but we give here 
a simpler direct proof which encompasses a more efficient algorithm, where we need 
to solve a program in 2d variables instead of nd variables. So we have the following. 

Corollary 3.1 The intersection problem (Q|) over any totally unimodular matrix A, 
any integer right-hand side b, and any positive integer n, is polynomial time solvable. 

Proof. Consider the program 

min{|z| : y G {0,1,..., n — l} d , £ G {0, l} d , A(y + z) = nb} . (7) 

We solve it in polynomial time using linear programming. If it is infeasible then so 
is program (j6j) and we are done. So assume we obtain an optimal solution y, z. 

Applying Lemma [2.11 to x := y + z, we can find in polynomial time x 1 ,... ,x n 
feasible in (J6]) with x k = x = y + z. For 1 < j < d, if xj = ■■■= xf = 1 

then yj + z 3 = J2k=i = n an d Vj < n — 1 so z 3 — 1. Therefore x 1 A • • • A x n < z. 
We claim the x k are optimal for (EJ). Suppose indirectly x 1 ,... ,x n form a better 
solution. Then z := x 1 A • • • A x n and y := Y^k=i % k ~ z are feasible in (J7|). But then 

\z\ = |a: 1 A---Ax n | < |x 1 A---Ax n | < \z\ , 
which is a contradiction. So x 1 ,...,!" are optimal and the corollary follows. □ 

Finally, we point out two variants of problem (EJ) where even deciding if the 
optimal value in the problem is 0 is NP-complete, even for fixed n = 2. First, if for 
each k we have a system Ax k = b k with the same matrix A but its own right-hand 
side b k , then, this is hard even if A is the vertex-edge incidence matrix of a directed 
graph G. Indeed, for k — 1,2, let Sk and A be distinct vertices in G and let b k be the 
vector forcing x k G {0, l} d satisfying Ax k = b k to be a path from sy to t^. Then the 
optimal objective function value of problem (15|) is 0 if and only if for k — 1,2 there 
are paths from sy to tk which are edge-disjoint, which is a classical NP-complete 
problem to decide. Second, if A is not totally unimodular, then this is hard even if 
A has a single row. Indeed, let ai,..., be given positive integers with even sum, 
and consider the NP-complete problem of deciding if there is a partition, that is, 
I C {1,..., d} with ^Zia a i = Yji^i a i■ Let A := ■ ■ ■, a d ], let b := \ Yf, >= \ a u and 

let n 2. Then, if / provides a partition, then x 1 , x 2 defined by xj := 1 — xf := 1 if 
i G / and xj := 1 — xf := 0 if i ^ /, are feasible in (Jl]) with value 0, and if x l , x 2 are 
feasible in (pQ) with value 0, then / := supp(x 1 ) := {i : xj ^ 0} provides a partition. 
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